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Effect of magnetic boundary conditions on the dynamo threshold 
of von Karman swirling flows 
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Abstract. - We study the effect of different boundary conditions on the kinematic dynamo thresh- 
old of von Karman type swirling flows in a cylindrical geometry. Using an analytical test flow, we 
model different boundary conditions: insulating walls all over the flow, effect of sodium at rest on 
the cylinder side boundary, effect of sodium behind the impellers, effect of impellers or side wall 
made of a high magnetic permeability material. We find that using high magnetic permeability 
boundary conditions decreases the dynamo threshold, the minimum being achieved when they are 
implemented all over the flow. 



> 

m 

o 

oo 

O 



% 



Dynamo action, i.e., self-generation of magnetic field 
by the flow of an electrically conducting fluid, is at the 
origin of planetary, stellar and galactic fields [1]. Fluid 
dynamos have been observed only recently in laboratory 
experiments in Karlsruhe [2] and Riga [3] by geometri- 
cally constraining the flow lines in order to mimic laminar 
flows that were known analytically for their dynamo effl- 
oiency [4]. More recently, the VKS experiment displayed 
self-generation in a less constrained geometry, e.g. a 
von Karman swirling flow generated between two counter- 
rotating impellers in a cylinder [5]. However, until now, 
dynamo action in the VKS geometry has been found only 
when the impellers are made of soft iron. It is thus of 
primary importance to understand how the dynamo prob- 
lem is modified by the presence of magnetic material at 
the flow boundaries. We address this problem here using 
a kinematic dynamo code in a cylindrical geometry. Two 
important approximations are made to simplify the study. 
First, an analytic test flow that mimics the geometry of 
the mean flow of the VKS experiment is considered. Sec- 
ond, the magnetic boundary conditions are taken in the 
limit of inflnite magnetic permeability of the boundaries 
compared to the one of the fluid. This seems a reasonable 
approximation for soft iron compared to liquid sodium. 
Our main result is that the critical magnetic Reynolds 
number, Rnic, for dynamo generation is significantly de- 
creased with boundaries of high magnetic permeability all 



over the flow. 

The VKS experimental set-up is sketched in Figure 1. A 
turbulent von Karman flow of liquid sodium is generated 
by two counter-rotating impellers (rotation frequencies Fi 
and F2). The impellers are made of iron disks of radius 
154 mm, fitted with 8 iron blades of height 41.2 mm, and 
are placed 371 mm apart in an inner cylinder of radius 
206 mm and length 524 mm. It is surrounded by sodium 
at rest in another concentric cylindrical vessel, 578 mm in 
inner diameter. This has been shown to decrease the dy- 
namo threshold in kinematic computations based on the 
mean flow velocity [6] . When the impellers are operated at 
equal and opposite rotation rates F, a statistically station- 
ary magnetic field is generated above a magnetic Reynolds 
number Rm ~ 30 [5]. The large scale field involves an az- 
imuthal component and a poloidal one which is dominated 
by an axial dipole. This geometry has been understood 
with a simple a — lu dynamo model [7] by taking into ac- 
count the helical nature of the flow that is ejected by the 
centrifugal force close to each impeller between succes- 
sive blades. Relying on the mean flow alone to compute 
the kinematic dynamo, smoothes out these non axisym- 
metric velocity fluctuations and thus cannot generate an 
axisymmetric field according to Cowling theorem. A non 
axisymmetric field is obtained, dominated by an equato- 
rial dipole [6,8]. 

When the disks are counter-rotating at the same fre- 
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Fig. 1: Sketch of the VKS experiment [5]. The inner and outer 
cylinders are made of copper (in gray). The dimension are 
given in millimeter (left) and normalized by the radius of the 
inner cylinder (right). 



qucncy, the structure of the mean flow (averaged in time) 
has the following characteristics: the fluid is ejected ra- 
dially from the disks by the centrifugal force and loops 
back towards the axis in the mid-plane between the im- 
pellers. A strong differential rotation is superimposed on 
this poloidal flow, which generates a high shear in the 
mid-plane. We approximate the experimental configura- 
tion with impellers of radius R (i.e. extending up to the 
inner cylinder boundary). We use cylindrical coordinates 
(s, (j), z). In this cylindrical domain ([0, 1] x [0, 27r] x [—1, 1]), 
the flow is well described by the following analytical ex- 
pression 

u = V X {tp e^) + su e^ , (la) 

with angular velocity lo and recirculation -0 respectively 
given by 
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This flow has been shown to generate a similar kinematic 
dynamo as the one computed using the experimentally 
measured mean flow [9]. In the above expression, e is a 
parameter controlling the ratio between the poloidal and 
toroidal components of the flow. We take e ~ 0.7259 as in 
previous numerical studies using this flow, in order to min- 
imize the critical magnetic Reynolds number for dynamo 
threshold [9,10]. 

Although we are aware that the experimentally observed 
dynamo cannot be captured with a kinematic calculation 
using the mean flow alone, we use this simple model here 
in order to study the effect of the magnetic boundary con- 
ditions on the dynamo threshold. We hope that the qual- 
itative behaviors will be unchanged with other dynamo 
modes in the presence of a turbulent flow. 



We perform direct numerical simulations of kinematic 
dynamos, solving the induction equation governing the 
evolution of the solenoidal magnetic field B 



-— ^ RmV X {u X B) + AB , 



(2) 



written in dimensionless form, using the diffusive 
timescale. The magnetic Reynolds number Rm is defined 
as Rm — fJ'OO'RUmax, where R is the cylindrical radius of 
the domain of the flow defined by (la-c) and Umax is the 
peak velocity of the flow. 

Equation (2) with the flow given by (la-c) is solved us- 
ing a finite volunre code adapted from [11]. To circumvent 
the severe CFL restriction induced by cylindrical coordi- 
nates, we ensure numerical stability near the axis using a 
low pass Fourier filter in the 0-direction. Also a centered 
second order scheme has been preferred fo to an up-wind 
scheme to discretize the inductive term, as resistive effects 
are here important enough to regularise the solution. As 
in [11], we ensure that V-B — is exactly satisfied using a 
constraint transport algorithm. The finite volume solver is 
fully three-dimensional. We have not used the decoupling 
of Fourier modes in the 0-direction. The initial magnetic 
field obviously needs to satisfy the divergence-free con- 
straint as well as the boundary conditions. In practice, an 
arbitrary divergent free field is initialized away from all 
boundaries. 

We investigate several types of magnetic boundary con- 
ditions. The classical approach is to use insulating bound- 
aries, matching the internal nragnetic field to the vacuum 
magnetic potential. The continuity of the magnetic field 
results in a non-local set of boundary conditions, which 
can be expressed via a "Neumann to Dirichlet" operator. 
We rely here on such an approach, using the boundary 
element formalism, as introduced in [12]. 

We investigate the effects introduced by using ferro- 
magnetic boundaries. This boundary condition can be 
expressed in a local form in the limit of infinite perme- 
ability. Jump conditions at a boundary between media of 
different magnetic permeability are well established 
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B X nl2 = — B X nji , 
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(3) 



(4) 



where subscripts 1 and 2 denote the two different regions 
and n is normal to the boundary. Ferromagnetic disks 
yield /i ^ /io, one can therefore reasonably approximate 
these jump relations, by using boundary conditions on the 
fluid side of the form 



B X n = 0. 



(5) 



This set of boundary conditions trivially implies that nor- 
mal currents vanish. 
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We will consider the effect of iron disks in two differ- 
ent ways. The first, and probably simpler, approach is to 
assume that the field is normal to the disks, namely 



B xe. 



0. 



(6) 



on the top and bottom of the cylinder. This boundary 
condition is well known in magnetohydrodynamics, in par- 
ticular in the astrophysical community, and is sufficient to 
close the system of equations we investigate. Another, and 
more subtle approach, allows to take into account the ef- 
fect of the blades on the disks. Assuming that the field 
becomes normal to radial blades as it approaches the end 
of the cylinder yields 



-B X Grf 



0. 



(7) 



This is an extremely simplified approach which does not 
take into account each individual blade, but accounts for 
their average effect on the large scale field. This set of 
boundary conditions is far less common, and deserves some 
care to ensure it provides the required constraints on the 
field. The solenoidal nature of magnetic field (V • S = 0), 
together with the fact that B^ is identically zero on the 
boundary, then imply 



dBd, 



BE, 

dz 



This yields the solvability condition 



27r 



OB, 
dz 



d<f> = 0. 



(8) 



(9) 



As the flow we consider here is axisymmetric, modes in 
the (/)-direction decouple, and we know from Cowling's 
theorem that the magnetic eigenmode cannot be axisym- 
metric. The solvability condition is therefore obviously 
satisfied (this would not be the case for a non-linear sim- 
ulation). 

For non axisymmetric modes, one can write 



B^dcl, = 0, 



(10) 



which together with (8) determines the B^ field com- 
pletely. The numerical implementation is two fold, we 
first compute 



B'M) 



dB^ 

dz 



(11) 



up to an arbitrary constant. Then we correct this function 
setting by the constant to meet (10). 

This implementation is mathematically consistent and 
maintains the divergence free property of the magnetic 
field. We shall however stress again that this very ideal- 
ized condition can only be used here because the modes 
decouple in the azimuthal direction and the axisymmetric 
mode cannot be unstable. 



The convergence of the numerical implementation has 
been carefully validated comparing simulations at different 
resolutions. We report results obtained with a resolution 
of 200 X 200 X 256 points. When an additional domain 
of sodium at rest is included, we use the same resolution 
in the inner flow domain defined by (la-c), which leads 
for the full domain to 240 points in the radial direction 
and 250 points in the z-direction to maintain a uniform 
accuracy. 




Fig. 2: The magnetic eigenmode obtained with potential 
boundary conditions. The rotation axis is vertical. An iso- 
surface of the magnetic energy (25% of the maximum value) is 
represented on the left. Magnetic field lines are plotted on the 
right. 

As in previous studies, all simulations yield magnetic 
eigenmodes with an m = 1 azimuthal symmetry. The 
structure of this eigenmode essentially corresponds to an 
equatorial dipole. This mode is represented on Figure 2. 

Previous numerical studies [13,14] compared the thresh- 
old values in configurations including and excluding the 
effect of fluid behind the disks. We reproduce here a simi- 
lar behaviour of the threshold: using insulating boundary 
conditions directly on the disks yields a threshold value 
Rrric = 63. This last value is in good agreement with [14]. 
Adding a layer of thickness 0.25 of fluid at rest between 
each disk and insulating boundaries increases the dynamo 
onset up to Rrric = 72. This is illustrated on Figure 3. 
Simulations, not discussed here, show that if the sodium 
behind the disks is not at rest, the threshold is even further 
increased. 

When we consider the case of ferromagnetic disks, the 
threshold for dynamo action is Rrric = 60.5 for a magnetic 
field normal to the disks and Rnic = 58 for a magnetic 
field in the i/i-direction. These two thresholds are very 
close and lower than in the case of a potential field (see 
also Figure 3). We thus observe that the high permeability 
boundary condition on the disks appears to do more than 
just screening the inhibition of electromagnetic induction 
due to the flow behind the disks. 

Implementation of a cylindrical layer of sodium at rest 
on the side can easily be included in the simulations. Since 
all electrical quantities are continuous, and in particular 
there is no jump in the conductivity, this does not involve 
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Table 1: Dynamo thresholds for different types of boundary 
conditions on the disks and on the sides of the cylinder. _B„ and 
B^i, denote ferromagnetic boundaries (respectively normal to 
the disks and normal to the blades), — V$ denotes the vacuum 
condition and Na indicates the presence of a layer of sodium 
at rest between the flow and the boundary. For the disks, cases 
Na & B^, Na & Bn or Na & — V<1? correspond to the presence of 
a layer of sodium at rest, outside which the relevant boundary 
condition is applied. The sodium lies behind the disks, and the 
boundary conditions are implemented on the lids. 



Fig. 3: Growth rate of the magnetic field as a function of Rm 
for different types of boundary conditions. We can see that 
replacing sodium at rest behind the disks by vacuum leads to 
a reduction of the dynamo threshold. Moreover, ferromagnetic 
conditions lead to additional reduction of threshold (about 10% 
for azimuthal conditions). Note that using high magnetic per- 
meability boundary conditions on the disks but also on the side 
yields the best configuration, with an onset at Rmc ~ 41.5. 



any extra boundary condition. In some sense, the model 
assumes a continuous fluid whose velocity field goes to 
zero outside a given domain. Table 1 reports critical pa- 
rameters calculated for different combinations of boundary 
conditions on the side and on the disks. The sodium at 
rest leads to an important reduction of dynamo thresholds 
when it is added on the side. It is known [10, 13, 14] that 
this effect increases with the width of the layer, but sat- 
urates relatively rapidly. In all the simulations presented, 
the width of the layer of sodium at rest on the side is set 
to 0.4 in order to study a configuration close to the VKS 
experiment. This positive impact on the onset appears to 
be independent of the boundary conditions on the disk. 

Finally, observing the strong effect produced by ferro- 
magnetic disks, it is tempting to investigate the effect of 
high permeability material for the entire vessel. We thus 
perform simulations for which the magnetic field is nor- 
mal to all boundaries (see Figure 3). This configuration 
appears to be the most efficient and it leads to a criti- 
cal magnetic Reynolds number Rnic = 41.5. It is clear 
from Table 1 that in the presence of sodium at rest, ferro- 
magnetic boundaries still yield a strong reduction of the 
onset, despite the fact that the boundary is now remote 
from the flow domain. No signiflcant modification of the 
global magnetic structure can be observed in our simula- 
tions when we compare ferromagnetic and vacuum bound- 



ary conditions. The magnetic eigenmodes in the bulk are 
similar with both sets of boundary conditions. 

The effect of the electrical conductivity of the bound- 
aries on the dynamo threshold has been studied since a 
long time [15]. It has been shown that the addition of a 
layer of electrically conducting fluid at rest around the flow 
can either decrease or increase the threshold [16]. This has 
also been observed in the case of von Karman flows, de- 
pending on the location of the layer at rest [13, 14]. Thus, 
general rules for the dependence of the dynamo threshold 
on the electrical conductivity of the boundary do not seem 
to exist. Such may not be the case with the magnetic per- 
meability of the bounding domain. Previous studies on 
the Ponomarenko dynamo with high magnetic permeabil- 
ity boundaries have displayed a decrease of the thresh- 
old [17, 18]. A similar improvement has been observed 
in the case of convectively driven dynamos [19]. In the 
present study, we flnd that boundaries with a high mag- 
netic permeability always decrease the dynamo threshold 
whatever their location. 

Two important aspects of the VKS experiment are not 
taken into account in the present study. First, the magne- 
tization of iron can lead to an additional amplifying factor 
of the dynamo as discussed in [7]. However, the coercitive 
fleld of pure iron being much smaller than the flelds gener- 
ated by the dynamo, the iron disks do not impose any per- 
manent magnetization. Reversals of the generated mag- 
netic field are indeed observed [20]. Second, the geometry 
of the magnetic field generated in the VKS experiment dif- 
fers from the one computed from the axisymmetric mean 
flow. Non axisymmetric velocity fluctuations generate a 
poloidal field with a dominant axial dipole, together with 
a strong azimuthal component. Thus the poloidal field is 
roughly normal to the disks, whereas the toroidal field is 
normal to the blades. It would therefore be of interest to 
check whether the dynamo threshold can be reached by 
using only ferromagnetic disks or blades. These two cases 
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would correspond respectively to our boundary conditions 
(6) and (7). 

Another interesting set of modifications suggested by 
this numerical study, would be to replace the copper side 
wall by a ferromagnetic one in the VKS experiment. If the 
experimentally realised dynamo mode behaves similarly 
to the simulations, we expect the threshold to be even 
lower in such configuration than that associated with only 
ferromagnetic disks and blades. 
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